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ABSTRACT 

"Smart"  structures  with  integrated  sensors,  actuators, 
and  control  electronics  are  of  importance  to  next-g‘-neration 
high-performance  structural  systems.  Piezoelectnc  materials 
possess  unique  electromechanical  properties,  the  direct  and 
converse  effects,  which  can  be  used  in  sensor  and  actuator 
applications.  In  this  study,  piezothermoelastic  charactenstics  of 
piezoelectric  shell  continua  are  studied  and  applications  of  the 
theory  to  active  structures  in  sensing  and  control  are  discussed. 
A  generic  piezothermoelastic  shell  theory  for  thin  piezoelectric 
shells  IS  derived  usinf  the  linear  piezoelectric  theory  and 
Kirchhoff-Love  assumptions.  It  shows  that  the  dynamic 
equations,  in  three  principal  directions,  include  thermal  induced 
loads  as  well  as  conventional  electric  and  mechamcal  loads.  The 
electric  membrane  forces  and  moments  induced  by  the  converse 
effect  can  be  used  to  control  the  thermal  and  mechanical  loads. 
A  simplification  procedure,  based  on  Lame’s  parameters  and 
radii  of  curvatures,  is  proposed  and  applications  of  the  theory  to 
1)  a  piezoelectnc  cylindncal  shell  and  2)  a  piezoelectric  beam 
are  demonstrated. 


INTRODUCTION 

Development  of  "smart"  structures  with  integrated 
sensors,  actuators,  and  control  electronics  are  crucial  to 
next-generation  structural  systems.  New  sensor/actuator 
materials  are  investigated  and  new  technologies  are  developed  in 
recent  years.  Among  those  commonly  used  sensor/actuator 
materials  (e.g.,  piezoSectric  matenals.  shape-memory  alloys, 
electrorheological  fluids,  electrostnctive  matenals. 
magnetostnctive  materials,  etc.),  piezoelectric  matenals  possess 
unique  electromechanical  properties  (the  direct  and  converse 
piezoelectnc  effects)  which  can  be  respectively  used  in  sensor 
and  actuator  applications  (Tzou  k  Fukuda.  1991:  Tzou  k 
Anderson.  1992). 


General  theories  derived  from  a  generic  shell  continuum 
can  be  applied  to  a  broad  class  shell  and  non-shell  structures 
(Soedel,  1981).  Chau  (1986)  proposed  a  vanational  formulation 
to  describe  the  electromechamcal  equilibnum  of  completely 
anisotropic  piezoelectric  sheUs.  Rogacheva 

(1982,1984a.l984b.l986)  studied  state  equations  and  boundary 
conditions  of  piezoelectnc  shells  polarized  along  coordinate 
directions.  Senik  and  Kudriavtsev  (1980)  formulated  the 
equations  of  motion  for  piezoelearic  shells  transversely 

polarized  Dokmeci  (1978)  derived  a  theory  for  coated 
thermopiezoelectric  laminae.  Tzou  and  Gadre  (1989)  proposed  a 
generic  theory  for  multi-layered  piezoelectric  shell  actuators 
based  on  equivalent  induced  strains.  Tzou  (1991)  denved  a 
general  distributed  sensing  and  control  theory  for  a  genenc  shell 
continuum  using  piezoelectric  thin  layers.  A  thin  piezoelectric 
solid  finite  element  with  three  internal  degrees  of  freedom  was 
formulated  and  applied  to  distributed  sensing  and  control  f 
continua  (Tzou  k  Tseng,  1991).  Tzou  and  Zhong  (1990)  denvir: 
a  piezoelastic  vibration  theory  for  a  hexagonal  symmetrical 
piezoelectric  thick  shell  with  three  effective  prinapal  axes;  and 
this  theory  was  applied  to  distributed  shell  convolving  sensors 
(Tzou  k  Zhong.  1991a)  and  active  structural  control  (Tzou  s: 
Zhong,  1991b)  In  this  study,  the  piezoelastic  shell  vibration 
theory  is  extended  to  include  thermal  induced  effects  due  to 
temperature  vanations.  Piezothermoelastic  behaviors  of 
piezoelectnc  shell  continua  are  investigated. 

Based  on  the  linear  piezoelectnc  theory  and 
Kirchhoff-Love  assumptions,  a  genenc  piezothermoelastic  shell 
theory  for  thin  piezoelectnc  shells  is  denved  first  .\ 
simplification  procedure,  based  on  Lame’s  parameters  and  radii 
of  curvatures,  is  proposed  and  applications  of  the  theory  to  a 
number  of  piezoelectnc  continua  (a  piezoelectric  cylindrical  shell 
and  a  piezoelectric  beam)  are  demonstrated.  Thermal  effects  to 
sensing  and  control  are  discussed. 
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It  is  assumed  that  a  genenc  piezoelectnc  shell  continuum 
IS  defined  in  a  curvilinear  tri-orthogonal  coordinate  system  m 


which  th«  ai  and  a:  define  the  neutraJ  surface  and  ai  defines  the 
normal,  Fieure  1  Since  the  shell  is  thm  the  electric  field  F3  is 
considered  across  the  shell  thickness  and  the  external  electric 
charge  Qj  is  on  the  top  and  bottom  surfaces  only  In  this 
section,  assumptions  and  constitutive  equations  are  defined, 
(Note  that  this  shell  is  generic,  which  can  be  simplified  to  a 
broad  class  of  shell  and  non— shell  geometries.  Examples  are 
demonstrated  in  case  studies,) 


f 


Fig,l  A  piesoelectric  shell  continuum. 

The  constitutive  equation  of  piezothermoelasticity  is 
defined  as 


It  IS  assumed  that  the  piezoelectric  shell  is  thin  as 
compared  with  the  other  two  in-plane  dimensions.  The 
transverse  shear  deformations  and  rotary  inertias  are  negligible 
Thus,  the  displacement  (r,  1  =  1,2)  of  any  given  point  in  the 

shell  continuum  can  be  represented  as  a  summation  of  the 
component  due  to  contraction/expansion  of  the  neutral  surface 
and  the  component  due  to  bending: 

Oi.rti.oii)  =  U|{ai,a2)  +  03^1(01,05), 

1  =  1,2,3.  (6) 


where  3^  denotes  the  bending  angle  and  J3  =  0  03  defines  tne 

distance  measured  from  the  neutral  surface.  Based  on 
Kirchhoff-Love  assumptions,  the  transverse  shear  strains  S.3 
and  S23  are  negligible,  i.e  .  S13  =  0  and  S23  =  0.  Thus,  the  two 
bending  angles  can  be  derived  as; 


3  —  ^  I 

3  _  U2 


1  <?U3 

x;  3o, ’ 

1  I?tl3 

Ts  3o2 


(') 

(8) 


Note  that  the  transverse  displacement  U3  is  independent  of 
thickness,  i.e.,  L''3  =  U3(0|.u2)  auu  the  transverse  strain  S33  can 
thus  be  neglected,  except  where  a  concentrated  load  is  applied. 
The  mechaiucal  strains  of  the  thin  shell  consist  of  an  in— plane 

membrane  strain  component  5°^  and  an  out-of— plane  bending 

component  k^j. 


where  {T}  is  a  stress  vector;  [c|  is  the  elastic  moduli  matrix;  [e] 

IS  the  piezoelectric  constant;  {)>}  =  [sj  is  the  elastic 

compliance  matrix;  {o’}  is  the  coefficient  of  thermal  expansion; 
{D}  is  the  electric  displacement  vector;  {S}  is  the  mechamcal 
strain  vector;  [<)  is  the  dielectric  constant  matnx:  {E}  is  the 
electric  field  vector;  {p}  is  the  pyroelectric  constant:  and  ,ltp  is 
the  temperature  change.  It  is  assumed  that  the 

piezothermoelastic  behaviors  are  instantly  balanced  in 
mechamcal,  electric,  and  thermal  fields  and  a  quasi-static 
approximation  can  be  applied.  For  a  piezoelectric  shell  with  a 
hexagonal  symmetncal  structure  (class  Cg^=  6mm),  the  elastic 

moduli  [c]  matrix  is  defined  by 


C11C12C13  0  0  0 

C|2C||C,3  0  11  0 

c  1 3  c  1 3  c  3  3  0  0  0 

0  0  0  c<4  0  0 

0  0  0  0  c«4  0 

.0  0  0  0  0  Cft, 


(3) 


where  Cn  —  (Y/1— u*),  cu  —  (Yji/1— 4J),  cjj - }-(C||-C|2) 

=  [Y/2( l-<-u)J  ,  where  u  is  Poisson's  ratio  and  Y  is  Young’s 
modulus.  (^Note  that  C|3,  C33,  C44  are  neglected  for  thin 
piezoelectric  shells  with  ineffective  in-plane  shear  constants.) 
Piezoelectric  constant  [ej  and  dielectric  constant  [<j  matrices  are 
defined  by 


r  0  0  0  0  0  0 1 

;e  ;  =  j  0  0  0  e;5e.,u  1 

^  1.  p!3  ei3  Cll  0  II  0  j 


1 


'  >  ,  I  0  o' 

!  0  Cil  I) 

L  0  0  4,3  J 


l41 


(5) 


Sii  =  S'))  +  Orikii  ,  (9~i) 

S22  =  S?:  +  oikw  ,  (9-t)) 

S|2  =  S'}:  +  oiki: .  {9-c) 


The  membrane  and  bending  strains,  S?j  and  k^,  are  defined  as 
follows: 


5}.= 
S?:  = 
5}:  = 

kii  = 

k::  = 

kt:  = 


1  ^1 1  u  2  3A|  ^ 

dO]  A I  A:  (fo:  H| ' 


I  C7U I  A I  A:  doi 

1  (7U2  ^  u  I  dXi 


SI' 


^2  ^^2  A  1 A2  C^oTj 
A2  d  ru2l  .  A,  5  Tu,] 

3  2  d\  1 


\^d3,  ^ 
1  ddz 


O  2  OAi 

A I A2  doj ' 
3 1  dA'i 
■A  I  .A 2  da  1  ' 


A:  ^.1^21  +  Aa  5  [/?,! 
"XJ  <7ai|.XJ  X  do2[Xj 


(10-a) 

(10-b) 

(10^) 

(11-a) 

(11-b) 

(ll-c) 


where  3s  are  defined  in  Eqs.(7)  and  (8).  Note  that  there  is  no 
shear  strain  on  the  03  face  such  that  there  is  no  induced  electnc 
field  in  the  Oi  and  the  03  directions.  Considering  the 
piezothermoelastic  cinstitutive  equations  and  the  stress-strain 

relations  of  thin  shells,  one  can  define  the  mechanical  stress  T^^ 
induced  by  the  mechamcal  strains,  the  electnc  displacement  Sj 
induced  by  strains,  and  the  stress  E|  induced  by  electnc  fields. 


Tvv  =  Ci3S3-.  -  C12S22  .  (l2-al 

T22  1=  Ci255i  —  C11S23  ,  (12— b) 

Tjj  =  0  .  (12-cl 

T|2  =  C4jS|2  (12-di 
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Ti3  =  ^23  =  0  . 

(12-e'i 

S,  =  S2  =  0 , 

(13-a) 

S3  =  eiiSii  +  631822 , 

(T3-bl 

El  =  Es  =  631E3 , 

(U-a) 

E3  =  633E3  , 

(14-b) 

E4  =  Eo  =  Ej  =  0  . 

(ll-c) 

where  the  first  term  is  contributed  by  the  converse  effect  :r.e 
second  term  by  the  pyroelectric  effect  I  temperature),  the  tn;rd 
term  by  the  elastic  strains  via  the  direct  effect  The  electric 
bending  moments  are 


M!i 


eiiEjajdas 

J  03 


h^  £ilj 

12  (33 


k 


11 


k 32 )  . 


These  terms  will  be  used  in  conjunctions  with  the  energy- 
expressions  and  the  vanational  equations 


FORCES  AND  MOMENTS 


M52  =  eiiEsOidas 
■I  03 


63  1 


f33 


-(kii-rk::) 


In  this  section,  all  forces  and  moments  introduced  by 
mechanical,  electric,  and  thermal  effects  are  defined.  These 
foree  and  moment  components  will  be  used  in  Hamilton's 
equation  when  deriving  the  shell  piezothermoelastic  equations, 
The  mechanical  membrane  forces  are 

N7i  =  f  Tiido3  =  A’(S'i i+iiS Jr)  ,  (15— a) 

J  03 

NIj  =  [  T22d03  =  A(S32+ixS],)  ,  (15-b) 

■'  03 

N72=[  T,2da3=  =  (13-<) 

■'03 


where  K  =  \'h/(l— n*)  is  the  membrane  stiffness  and  Njj  is  the 
total  force  acting  on  the  ith  face  in  the  jth  direction  due  to 
mechamcal  effects.  The  mechanical  bending  moments  are 


M*i  =  [  TnOidoi  =  Dfkii-riik::)  ,  (13— a) 

■*03 

M?2  =  [  Tjja3dQ3  =  Z?(k22-+-Mi<n)  .  (16— b) 

■'  a3_ 

M72=f  T,203da3  =  Mil  =  ,  (16-<) 

■’03 

M73=M73  =  0,  (16^) 


Vh  3  ® 

where  D  =  ^2(  1-J7]  bending  stiffness  and  M  jj  is  the  total 
bending  moment  on  the  ith  face  in  the  jth  direction  due  to  the 

mechanical  effects.  The  mechanical  transverse  shear  forces  Qis 
are 

Q73  =  [  Tiidoi,  (1”~^) 

Ja3 

Q?3  =  [  Tsidoi .  (!'— b) 

Jq3 

Using  Eq.(2),  one  can  derive  the  electric  membrane  forces: 


Mfj  =  .\I'[3  =  M53  =  0  , 


where  M7j  is  the  total  electric  bending  moment  on  the  ith  face  :n 
the  jth  direction  due  to  the  converse  piezoeiectnc  effect  The 

electric  transverse  shear  forces  are  Q73  =  QJs  =  0.  There  is  no 
shear  forces  in  the  03  direction  due  to  the  electric  effects  The 

t 

thermal  membrane  forces  Nji  are  defined  by 

—  [  AiAtpd(i3  —  hAiAtp  ,  (20— at 

NJj  =  j  AjAtpdai  =  hAiAtp  .  (20— bi 

■’as 

The  thermal  bending  moments  are 

Mil  =  [  AiAtpOidai  =  0  ,  (21-ai 

■’as 

Mis  =  [  AsAtpOidtts  =  0  .  (21-bi 

■'as 

It  is  noted  that  the  piezoelectric  continua  experience  only 
in-plane  thermal  expansion/contraction  and  no  bending 
moments  in  a  uniformly  distributed  thermal  field. 


HAMILTON’S  PRINCIPLE  AND  PIEZOTHERMOELASTIC 
EQUATIONS 

In  this  section,  piezothermoelastic  equations  in  three 
principal  directions  vnll  be  derived  using  Hamilton’s  principle 
and  the  vanational  procedures.  Hamilton’s  principle  gives 
(Tzou  &  Zhong,  1990) 


4^‘f  ;^-pUjU2-H(Skj,Ej 

to"^  \  “ 

-)-  [  [  '  fj  5Uj  -  ^^(p)dSdt  =  0  , 
to-'  S 


dVdt 


N7i  =  [  esiEsdas 

■’“s 

=  -|ilh(?3-^hp3Atp-^S1,-  Ssslh,  (l8-a) 

«33  <33  '  «33 

•Vis  =  f  e3iE3da3 
•'03 

=  -^h(?3  hpjAtp  -  j2^S'J,-rS;:)h  ,  (l8-b) 

*33  ^33  ^33 

N72  =  0.  (l8-c) 


where  H  is  the  electric  enthalpy;  p  is  the  mass  density;  Uj  is  the 

displacement:  E,  is  the  electric  field;  ti  is  the  surface  traction  ;n 
the  aj  direction:  Syj  is  the  strain  on  the  kth  face  and  in  the  jth 

direction:  ^  is  the  surface  charge;  and  >p  is  the  electnc  potential 
The  electric  fields  in  the  curvilinear  coordinate  system  are:  1 1  E  ., 

=  -A,(l-l-a,)R,)l^’2)  = 

=  - ,  where  Ri  and  Rs  are  the  radii  of  curvatures,  and  A; 

and  As  are  Lame’s  parameters.  These  define  the  electnc  fields 
as  the  gradient  of  the  electric  potential.  The  electnc  enthalpy  is 
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defined  as  H  =  j[{S}‘{T}  -  {E}‘{D}]  (Tzou  ic  Zhong,  1990). 

Using  the  piezothermoelastic  constitutive  equations,  one  can 
denve 


H  =  MS}‘[cl{S}-{E}‘[el{S}-j{E}'[q{E} 

-{S}‘{A}Atp-{E}‘{p}Atp.  (23) 

Substituting  the  strain— stress  expressions  into  the  electric 
enthalpy  gives 

H  =  2(T|iSii-^T;2Sj;+T];S  12)  —  63  .(S  1 1-1-S22)  E  3 
-  Kf33E3^)  —  (AiS!i+A2S22+A3S33-rP3E3)Atp  (24) 

Substituting  the  electric  enthalpy  and  all  other  energy 
expressions  into  Hamilton’s  equation  and  collecting  all  like 
terms  in  the  variational  equation,  one  can  denve  the 
piezothermoelastic  and  vibration  equations  in  three  principal 
directions. 


(25-a) 


-  (N?r-N5r-NW^^  =  nhA.A,'’" 


at‘ 


da  2 
(26-a) 


(27-a) 


where  h  is  the  thickness  of  piezoelectric  shell.  The  superscripts 
m,  e,  and  t  respectively  denote  mechanical,  electric,  and  thermal 

components.  Q“3  and  Q'Js  in  Eqs.(25)— (27)  are  defined  by 


II  ,  ^fMTiA.i 


Q,.A,A,-aM.TlAJ^at(M!,-M1,lA.i 


(28) 


(29) 


Note  that  the  the  equation  of  motions  include  the  mechanical 
forces/moments  electric  forces/moments  (N?yM^), 

and  thermal  induced  forces  (N‘  )  As  discussed  previously. 

thermal  expansions/contractions  in  three  principal  directions  are 
considered.  It  is  observed  that  the  umform  temperature 
variation  does  not  contnbute  any  thermal  moments,  which  will 
not  be  the  case  in  non-uniform  temperature  variations.  These 
system  equations  can  be  solved,  with  the  appropriate  boundarv 
conditions  and  external  excitations  (mechanical,  electric,  and/or 
thermal),  to  describe  the  exact  piezothermoelastic  behaviors  of 
the  piezoelectric  shell  In  structural  c.introl  applications,  the 
electric  related  components  can  be  used  as  control  forces  and 
moments  to  alter  system  characteristics  iTzou  Zhong.  1991b). 
Rearranging  Eqs  (25-27/  and  moving  all  electric  related  terms 
(control  terms)  to  the  right  (Tzou.  1991),  une  can  derive 
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Note  that  electric  related  terms  can  be  used  as  control  terms  to 
actively  change  the  shell  dynamics.  In  addition,  all  terms  with  a 
constant  either  1/Ri  or  l/Ri  vanish  if  the  radius  of  curvature  is 
infinite,  e.g.,  flat  plates  (Tzou,  1991).  The  charge  equation  of 
electrostatics  of  the  piezoelectric  shell  is  denved: 


i9f(e3iSii-(-e3|S22-)-t  3  3E3+P3AtD)A|A2l  _  g 
3q3 


(301 


which  implies  that  the  quantity  [e3i(Sii+S22)  +  fisEs  +  PsAtpj 
.AiA:  is  equal  to  a  constant  and  the  tUckness  variation  is  equal 
to  zero.  Note  that  this  equation  can  be  used  to  estimate  an 
electnc  output  as  functions  of  induced  mechanical  strains  and 
temperature  variation,  i.e..  E3  =  -  (l/t33)[e3i(Su+S22)+P3Atpi. 
in  an  open-circuit  condition  (Tzou,  Zhong,  1991a). 


BOUNDARY  CONDITIONS 

Boundary  conditions  are  directly  derived  from  the 
variational  equation  (.Appendix).  The  boundary  conditions  are 
defined  by  the  surface  traction  forces  and.  the  surface  charge 
(Note  that  other  types  of  boundary  forces  and  moments,  such  as 
spring  supported  boundaries,  fixed/hinged  boundaries,  etc.,  can 
also  be  accommodated.) 

Mechanical  Boundary  Conditions 

Mechamcal  boundary  conditions  defined  by  either 
force/moment  or  displacement/rotation  are  summarized  in 
Table  I  in  which  terms  with  a  superscript  *  denote  external 
boundary  components. 
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Table  1.  Mechanical  boundary  conditions. 


Fo  r  c  e  B.C  '  s 

Disp  B  C. 

>.m 
■  kk 

“  ^  kk  “  ^kk  ~  ■'’kk 

• 

^  =  “k 

\.ra 
k  k 

- 

\  =  -^k 

« 

« 

k3 

=  V,  0 
<3 

U]  =  U3 

Qkt 

lit  =  uj 

where  k  =  1.2.  the  subscript  t  denotes  the  taneential  direction 
(i.e..  t  =  2  if  k  =  1  and  vice  versai.  It  is  observed  that  the 
thermaJ  induced  membrane  force  only  occurs  ;n  the  principal 
direction.  The  mechamcai  shear  stress  resultants  are  defined  as 


\  13  -  Q  1 3  -  “ 

V.3  =  q'I'j - L  ,  (3l-a.b) 

Qi2  =  N’T:  ^  and 

Qj,  =  nTi  ^  1.32-a.bi 


.\gain.  there  is  no  electrically  induced  shear  components  because 
the  in-plane  twisting  effect  is  neglected.  Note  that  usually  ordy 
either  force  boundary  conditions  or  displacement  boundary 
conditions  are  selected  for  a  given  physical  boundary  condition. 

For  a  totally  fixed  edge  at  ai  =  0‘\  U  e  .  no  motion  allowed), 
the  boundary  conditions  are:  ui  =  0,  J,  =  0.  ui  =  0.  and  uj  =  0. 

For  a  totally  free  edge  at  03=02  ,  1  e  .  no  external  forces  and 
moments,  the  boundary  conditions  at  03  =  03  are:  N’Ts  —  N22 

=  0.  mT:  -  -Ms:  =  0.  V23  =  0.  and  T21  =  1'  In  the  case  where 
the  surface  traction  forces  t-^  are  defined,  the  boundary 


membrane  forces  are 

Nii=f  t,,do3, 

(33-a) 

■^03 

N22  =  t32d03  , 

{33-b) 

03 

N  12  =  ti3da3  , 

(33-c) 

03 

-N21  =  t2id03  , 

(33-d) 

J03 

where  N*j  is  the  total  force  on  the  ith  (ace  in  ihe  jth  direction 

due  to  the  surface  tractions.  The  induced  boundary  bending 

« 

moments  M  i  j  are 


VI 

M 

M 

M 


• 

1  1  = 

tn03d03  , 

(34-a) 

’  03 

27  = 

t2:a3dQ3  . 

i34-b) 

^  ai 

1  2  = 

t  ]:03da3  . 

(34-c) 

J  03 

• 

21  = 

tuaidoi 

(34-d) 

^  03 

.Vccordingiy,  the  boundary  transverse  shear  forces  Q  ,.  are 

Q 13  =  Q  3 1  ~  I  ^  iidos  .  .  .’5— a  1 

J  03 

•  *1' 

Q:3  =  Q32  =  t:3d03  .'ij-D. 

J  03 

where  Q  ,,  is  the  shear  force  cn  the  ith  face  in  the  jth  direction 


Electric  Boundary  Condition 

The  electric  boundary  condition  is  defined  as 

« 

e3lSii  -i-  e3|S22  -r  «33E3  +  P3Atn  "  (^3  =  ' 

It  is  observed  that  the  total  surface  charge  including  the 
mechamcai.  electric,  and  temperature  effects  is  equal  to  the 

external  surface  charge  (?3 

Note  that  the  piezothermoelastic  equations  for  the  thin 
shell  continuum  and  the  boundary  conditions  can  be  reduced  to 
conventional  elastic  shell  equations  by  neglecting  all  electric  and 
thermaJ  coupling  terms  (Soedel,  1981).  .Vgain.  transverse  shear 
deformation  and  rotatory  inertia  effects  were  not  considered 


PIEZOTHERMOELECTRICITY  OF  SIMPLIHED 
GEOMETRIES 

The  piezothermoelastic  theory  derived  above  is  for  a 
generic  piezoelectric  shell  continuum  exposed  to  mechaiucal. 
thermal,  and  electric  fields.  The  generic  shell  was  defined  in  a 
curvilinear  tri-orthogonal  coordinate  system  defined  by  a,,  o;. 
and  03  axes.  The  in— plane  two  axes  define  the  neutrjd  surface 
experiencing  only  membrane  effects.  Each  of  the  in-plane  a»s 
is  defined  by  its  radius  of  curvature,  e.g.,  Ri  for  oi  and  R:  for 
03.  In  addition,  there  are  two  Lame’s  parameters  (Ai  and  .Vi) 

2  2  2 

defined  by  a  fundamental  form:  (ds)  =  .4i‘  (doti 

+  .A3  (do3)  For  a  given  geometry,  Ri  and  Rs  can  usually  be 
directly  observed  from  the  coordinate  system  and  A,  and  .As  can 
be  derived  from  the  fundamental  form.  Substituting  the  four 
parameters  into  the  generic  shell  equation  and  simplifying  them 
accordingly,  one  can  derive  the  corresponding  piezothermoelastic 
equations  and  boundary  conditions  for  the  geometry  In  this 
section,  these  procedures  are  used  to  derive  the 
piezothermoelastic  equations  for  1)  a  piezoelectnc  cylindrical 
shell  and  2)  a  piezoelectnc  beam. 


Example— 1:  Piezoelectric  Cylindrical  Shell 

It  is  assumed  that  the  cylindncal  shell  is  defined  in  a 
cylindrical  coordinate  system  in  which  x  axis  (0|)  is  aligned  with 
the  height  and  its  radius  of  curvature  R|  =  a.  The  second  axis 
(03)  defines  the  circumferential  direction  which  has  a  radius  of 
curvature  Ri  =  S  Note  that  the  x  and  S  axes  constitute  the 
neutral  surface.  The  third  axis  03  is  normal  to  the  neutral 
surface.  Figure  2  illustrates  the  piezoelectric  cylinder  and  us 
coordinate  system.  Piezothermoelastic  effects  of  the  cylindrical 
shell  will  be  discussed. 
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Note  that  the  superscnpts  m.  e,  and  t  are  for  the  mechaiucai. 
electric,  and  thermal  effects  respectively  Substituting  the  force 

and  moment  terms  into  the  QTi  and  equations,  one  can 
denve 


r  r.  Oirh^’  'D(l-t-ai  ei'-h^l 

^ =  i “  ^  ■  777  Hi  ^  *  ^~ir- - 

r  D  ej,i  hJ  1  d3u, 
dddx  I  ~]n  t33  T2Tnj  'btdtl3  ■ 


I  4-1 1 


Qh  = 


r  D  e3i3  h3  4  <43113 
j_  S  3  6  33  12Ji3j  dS3 


.  '  D  e3i3h3  ’a3u,  D 

I ~  nr  632  ' nn 


Fig.2  A  piezoelectric  cylindrical  shell. 

The  fundamental  form  of  the  cylinder  is 

(ds)2  =  ll)3(dx)3 -(- S3fd«)2  .  '37) 

Thus,  .\i  =  1.  .^2  =  S.  Ri  =  I.  R:  =  S.  ol\  =  x,  a-i  =  0 
Substituting  these  parameters  into  the  member/ bending  strain 

expressions  in  Section-2,  one  can  denve  the  membrane  S?,  and 
bending  strains  kij  for  the  cylindrical  shell.  Thus,  the  total 
strains  are 


c  ^  ^^U3 

Su  =  ^-<"3^. 
1  ri3u«  1 


(38-ai 


c  -  1  [  *'  n  1  ^  1  f38-bl 

^53  =  nr  nrr  ■  (38  b) 

5  1  du,l  a,  JJui] 

Substituting  the  strains  into  the  stress  equations  and 
consequently  into  the  force/moraent  equations,  one  can  derive 
the  force/moment  for  the  cylindrical  shell: 


N--,=  [^-+-  uijl 


r  j  |-du, 


^'^3  =  K  [tt  [w  ^1  • 

Afi-uir^’ ,  1  ai,i 

^'3  =  ^3 — [w^'irw 


M"'  -  n\-  ^  I 

Mq  -  /O  r‘^‘'_53U3l  _  ^  ^ | 
'*32  -  ^  Jf7  1^  cfzrrj  4*^!  ■ 

Ml  =  -f^h  +  ^  U31! 

-  ^  *3  [pi-Ntp  . 

N«2  =  N'm 


(39-a) 

(39-b) 


(41-a) 

(41-b) 


\r  -  e3jjh3  <  d2ii3  ..  I  ^  Mil'll.  i  P-ji 

'‘‘■■  =  -777-n!:-'57r*-Tr  ,rr-7rr'7  ' 

Nl't;  =  Ml,  i42-bl 


Thus,  the  piezothermoelastic  equations  in  three  principal 
directions  for  the  piezoelectnc  cylindncal  shell  are  derived 

^(NL^NIcaH  .  SU  .  ,  h6, 


+  Q?3  =  ph®- 


_  rN?j_N5^.v5jl  =  Ph®^ .  (48 ) 

It  IS  observed  that  the  thermal  effects  only  contnbute  to  the 
membrane  forces.  Removing  the  electric  and  thermal  related 
terms,  one  can  simplify  the  system  equations  to  those 
corresponding  to  an  elastic  cylindrical  shell. 

Using  the  charge  boundary  condition,  one  can  define  the 
electric  field  strength  E3  at  the  location  03  above/ below  the 
neutral  surface  as  a  function  of  the  mechanical  strains, 
temperature  effect,  and  charge  effect. 


^3  =  ~  7^  I  ^  +  -i-[^  ^  “3]  ■  ‘‘3  ^ 

Dn  oifl  J  f  3  3  r  33 


The  electric  field  strength  is  contnbuted  by  the  direct 
piezoelectnc  effect  (the  first  term),  the  pyroelectnc  effect  (the 
second  term),  and  the  external  surface  charge  (the  third  termi 
as  defined  in  the  constitutive  equation.  Note  that  the  resulting 

voltage  is  V3  =  J*^^E3da3  in  an  open-circuit  condition.  The 

bending  components,  with  03  terms,  vanish  after  the 
integration.  It  is  also  observed  that  the  output  signal  has  a 
temperature  related  terra  induced  by  the  pyroelectric  effect  in 
sensor  applications  Note  that  it  is  assumed  that  the  external 
charge  is  zero  in  sensor  applications  (Tzou  k  Zhong,  1991a). 

Elxample-2'.  Piezoelectric  Beam 

.\  beam  is  a  special  case  of  an  open  nng  with  zero 
curvature.  5  =  i  In  this  case,  the  a,  axis  is  aligned  with  the 
longitudinal  direction  of  the  cantilever  beam,  1  e.,  q,  =  x.  The 
second  axis  is  in  the  width  direction,  ar  =  y  Figure  3  shows  the 
piezoelectnc  beam.  It  is  assumed  that  the  beam  only 
experiences  transverse  oscillations,  Q3  =  z  Goverrung  equation 
and  piezothermoelastic  behaviors  of  the  beam  are  discussed. 


M,  =  \:Al„ 

=  A.At,, 


(43-ai 

(43-bl 
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The  fundamental  form  of  the  beam  is 
(ds)3  =  (li-'(dx|3  -i-  (l)3(dy)3  . 


r 


Fig. 3  A  piezoelectnc  cantilever  beam. 

where  dx  and  dy  are  inCnuesimal  distances  in  the  x  and  y 
directions  respectively-  Thus,  .\i  =  i.  =  i.  Ri  =  i.  Ri  =  t 
Since  only  the  bending  oscillation  is  considered,  the  membrane 

strains  are  zeros,  i.e..  S'ii  =  0.  S::  =  <).  and  S';:  =  0  _  The 
bending  strain  at  the  03  location  is  defined  by  k.i  = - 

and  kjj  =  0,  k^  =  0.  The  total  strains  at  03  location  ate 
defined  as 

S„  =  -a3-^^.S2:  =  0,  S,;  =  '1  (51) 

.\gain,  the  beam  expenences  only  transverse  oscillation  The 
membrane  (longitudinal)  force  components  are  all  zeros,  i.e..  N” , 
=  0,  N'!:  =  0.  N”:  =  0-  The  resultant  moments  are 

M"i  =  Dfkii+iikjj)  =  -  D  ’  t,52~a) 

Ml2=  D(k:2+iiku)  =  (52~b) 

M72=  ^■2~-  l<i2  =  0  (52-c) 


Note  that  the  moment  M'J:  is  primarily  introduced  by  Poisson  s 
effect.  The  electnc  force  and  moment  resultants  due  to  the 
external  charge  and  temperature  are 


'■'ll  =  -  h  -^11—  Oz  -  h-^l-!_  0,  jitn  .  (oG-ai 

'  ‘33 

N«:  =  -h-^g3-h-^P3At^  (53-b) 

f33  <33  • 


NHi  = 

M?:  = 


h^  e3|i 
Tr-^-dxT- 
hJ  esi’  i92u3 


(54— a) 
'  54— b  I 


N(|  =  h  ,  (55- ai 

N):  =  h  A|Atp  l55— bi 


Substituting  the  system  parameters  and  force/moment 
resultants  into  the  original  shell  equation,  one  ran  derive  the 
transverse  piezothermoelastic  equation 


h^  e3|-’1  '^*113 
12  '33  j  chf* 


1 56) 


For  a  beam  with  a  rectangular  cross-section  1  width  b  and 
thickness  hi,  the  transverse  equation  of  motion  is 


M  - 


H'l 


'  13 


f!-  =  ohb 


,421 


Jt- 


where  I  =  — p-  Note  that  the  elasticity  pan  has  one  m  .re 

term  contributed  by  the  piezoelectricity  The  piezoeiectricitv 
■’ontributed  elasticity  is  very  small,  about  I'T  for  piezoelectric 
poiyvinylidene  fluoride  polymer  iTzcu  k  Zhong.  ;jy;i. 
However,  the  temperature  has  no  contribution  to  the  transverse 
oscillation  because  the  thermal  forces  are  primarily  n  •  r.e 
neutral  surface,  neutral  axis  in  this  case  This  pyroelectric  effect 
will  contribute  to  the  longitudinal  oscillation. 


The  electnc  field  strength  at  the  location  03  above 'bef.w 
the  neutral  axis  is  defined  by  the  external  charge,  temperature 
induced  pyroelectric  effect,  and  bending  strain 


Ej  =  -^ - ^P3-5tp 

f  33  tl3 


^31  .  '3‘IJ3! 


L‘33 


i'  'I 


However,  the  resultant  open-circuit  voltage  V3 
contnbuted  by  the  pyroelectnc  effect  I  C?3  =  0 


is.  in  fact. 

.h/2 

and  I 
'-h/2 


niy 

’3 


SUMMARY  AND  CONCLUSIONS 

■A  linear  piezothermoelastic  theory  of  piezoelectnc  shell 
continua  was  proposed  and  piezothermoelastic  phenomena  were 
evaluated.  It  was  assumed  that  the  electnc.  thermal,  and  elastic 
fields  are  instantaneously  balanced  and  a  quasi-static  >:  ndition 
IS  used  in  the  piezothermoelastic  constitutive  equations  A 
genenc  theory  for  a  piezoelectnc  thin  shell  continuum  .as 
denved  using  Hamilton’s  pnnciple  and  Kirchhoif  i.ove 
assumptions  The  govermng  equations  show  close  coupling 
effects  among  electnc.  thermal,  and  elastic  fields.  Both 
mechamcal  and  electric  effects  contnbute  to  the  resultant 
forces/raoments  for  the  shell  continuum.  However,  it  was 
observed  that  the  thermal  effect  only  contnbutes  to  the 
membrane  force  resultants,  not  the  bending  resultants  due  to  a 
uniform  temperature  assumption  Thermal  induced  bending 
could  appear  if  there  is  a  non-umform  temperature  distnbution. 
Note  that  the  electnc  fcrce/moment  resultants  in  the 
piezothermoelastic  equations  can  be  used  to  control  the  shell 
continuum. 

The  denved  piezothermoelastic  equations  are  genenc. 
which  can  be  simplified  to  a  vanety  of  piezoelectnc  continua  if 
two  radii  of  curvatures  and  two  Lame’s  parameters  are  defined 
This  simplification  was  demonstrated  in  three  examples;  D  a 
cylindncal  shell  and  2)  a  beam.  Detailed  piezothermoelastic 
phenomena  of  each  geometry  were  discussed  along  with  the 
denved  governing  equations.  The  same  procedure  can  be 
applied  to  a  variety  of  other  piezoelectnc  continua  and  so  as  the 
piezothermoela.snntv  evaluated.  .Note  that  the  theory  was 
denved  based  on  linear  assumptions  and  the  material 
nonlinearity  was  not  considered.  However,  these  material 
constants  (e  g  ,  piezoelectnc  constants,  elastic  constants,  etc  1 
could  vary  when  temperature  vanation  is  significant.  Thus, 
extending  the  present  theory  to  encompass  the  material 
nonlinearity  would  further  enhance  the  theoretical  development 
and  understand  more  about  the  complicated  behaviors  of 
piezoelectric  sensors/actuators  operating  in  non-ideal 
environments. 


211 


ACKNOWLEDGEMENT 


This  research  was  supported  by  a  i^rani  from  the 
National  Science  Foundation  (No  Ril-'*61067I  i  and  the 
Kentucky  EPSCoR  Program.  A  grant  from  the  Army  Research 
Office  1  DAAL03-91-G-<1065i,  Technical  Monitor  Dr  Gary  L 
Anderson,  is  also  gratefully  acknowledged  Contents  of  the 
information  do  not  necessanly  reflect  the  position  or  the  policy 
of  the  government,  and  nor  official  endorsement  should  be 
inferred  Invaluable  discussions  with  Professors  T  R  Tauchert 
and  G  E.  Blandford  are  also  gratefully  acknowledged 


REFERENCES 

Ddkmeci.  M  C  .  1978.  "Theory  'jf  Mbrations  of  Coated. 
Thermopiezoelectric  Laminae,"  J  Math  Phys  .  19i  1 1.  January 
Tlogacheva.  .N  N  .  1982,  'Equations  A  State  of 

Piezoceramic  Shells."  PM.M  C  S  S  R  .  i  pp  o77-^84 

Rogacheva,  .N  N  ,  1984.  'On  Stain-V<^nant  Type 
Conditions  in  the  Theory  of  Piezuelastic  Sheiis,"  PMM 
rSS  R  .  18(21,  pp  213-216’ 

Rogacheva,  N  .N  .  1984,  "On  Boundary  conditions  in  the 
Theory  of  Piezoceramic  Shells  Polarized  .^long  Coordinate 
Lines."  PMM  L’  S  S  R  .  47(2),  pp  220-226 

Rogacheva.  .N  N  1986.  "Classification  of  Free 
Piezoceramic  Shell  Vibrations,"  PMM  USSR.  501 1 1 

pp  106—1 1 1 

Senik.  N  .A  and  Kudriavtsev,  B  A  .  1980  "Equations  on 
the  Theory  of  Piezoceramic  Sheiis,  ’  In  Mechanics  of  a  soiid 
ieforrnable  body  and  related  analytical  vrohlems.  Moscow,  izd 
rnask.  Inst.  Chim.  .Afoshinostroenita.  U  S  S  R 

Tzou,  H  S  ,  1991,  "Distnbuted  .Modal  Identification  and 
Vibration  Control  of  Continua:  Theory  and  Applications  ' 
ASME  Journal  of  Dynamic  Systems.  Measurements,  and 
Control.  113(3).  pp  494— 199.  September  1991 

Tzou.  H  S  and  Anderson.  G  L  .  1992.  Intelligent 

Structural  Systems.  Kluwer  Academic  Pub  .  The  Netherlands 
t  To  appear) 

Tzou.  H  S  and  Fukuda,  T  ,  1991,  Piezoelectric  Smart 
Systems  .Applied  to  Robotics.  .Micro-Systems.  Identification,  and 
Control.  Workshop  .Notes.  IEEE  Robotics  and  .Automation 
Society,  1991  IEEE  International  Conference  on  Robotics  and 
.Automation.  Sacramento.  C.A.  .April  7-12.  1991 

Tzou.  H  S  .  and  Gadre.  M.  1989.  "Theoretical  .Analysis  of 
a  Multi-Layered  Thin  Shell  Coupled  with  Piezoelectric  Shell 
Actuators  for  Distnbuted  Vibration  Control."  Journal  of  Sound 
and  l  iiration.  132(3).  pp  433-450 


Tzou.  H  S  .  and  Tseng.  C  I  .  ;591,  "Distributed  .Mc-dai 
Identification  and  Vibration  Control  ■:{  C.mtinua  Piezoeiectr.c 
Finite  Element  Formulation  and  .Analysis.  '  ASME  Jouma.  • 
Dynamic  Systems  .Measurements  ana  ('onfoi.  IIJ.'I 

pp  500-505 

Tzou.  H  S  i;  Zhong.  J  P  .  liOO.  ElectromechaniraJ 
Dynamics  of  Piezoeiectnc  Shell  Distnbuted  Systems.  Part  : 
Theory  and  Part-2  Applications.  "  Robotics  Researcn- :  .o-i 
ASME-DSC-Vol  26.  pp  199-21  1  1990  ASME  AVinter  .Annuai 
.Meetings.  Dallas.  Texas.  Nov  25-!')  1990  tna 

"Electromechanics  and  Vibrations  A  Piezoeiectnc  'ihea 
Distnbuted  Systems."  .AS.ME  Journal  of  Dynamic  Systems 
.Measurements,  and  Control.  iTo  appear  i 

Tzou.  H  S  and  Zhong.  J  P  .  1991a.  "Sensor  .Mechanics  : 
Distributed  Shell  Convolving  Sensors  Applied  to  Flexible 
Rings."  Structural  Cibration  and  .Acoustics.  Edrs  Huang.  Tz  u 
et  al  .  .ASME-DE-VoI  34.  pp  67-74.  Symposium  on  Intellieer.t 
Structural  Systems.  1991  ASME  13th  Biennial  Conference  n 
Mechanical  A’lbration  anc  Noise.  Miami.  Florida.  Septemcer 
"22-25.  1991  .ASME  Journal  of  'Aibration  ana  .Acoustics  T  . 
appear ) 

Tzou.  H  S  and  Zhong,  J  P  ,  I991b.  "Control  A 

Piezoeiectnc  Cylindrical  Shells  via  Distnbuted  In-P!ane 
Membrane  Forces."  Controls  for  .Aerospace  Systems. 
DSC— Vol.35.  pp  15—20.  Distnbuted  Control  of  Flexible 
Structures,  1991  .ASME  W.A.M,  .Atlanta.  G.A.  December  1—6 
1991  '  PiThmiShl  WAM92 


212 


